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INTRODUCTION 

In  many  streams  and  rivers,  the  proportion  of  material  transported 
as  bedload  is  a  relatively  small  percentage  (<20%)  of  the  total  sediment 
load.  However,  as  the  bed  and  banks  of  a  "self-formed"  alluvial  channel 
are  composed  largely  of  this  material,  its  potential  for  movement 
ultimately  controls  the  stability  of  the  river.  Further  interest  in 
bedload  transport  phenomena  has  been  sparked  by  a  number  of  recent 
studies  which  demonstrate  high  spatial  and  temporal  variability  in 
transport  rates. 

Although  interest  in  the  topic  is  high,  complete  sets  of  bedload 
sample  data  are  relatively  scarce.  The  purpose  of  this  report  is  to 
present  the  results  of  sediment  discharge  measurements  made  at  Fall 
River,  Rocky  Mountain  National  Park  during  the  1983  summer  field  season. 

FIELD  SITE  AND  DATA  COLLECTION 

The  Lawn  Lake  flood  of  July  15,  1982  resulted  in  impressive  channel 
changes  along  Roaring  River  and  Fall  River,  Rocky  Mountain  National  Park 
(Figure  1).  Flood  waters  incised  deep  gulleys  into  the  glacial  deposits 
which  mantle  the  Roaring  River  valley.  Much  of  this  eroded  material  was 
then  deposited  as  a  coarse  debris  fan  in  the  broad  low-gradient  valley 
of  Fall  River  in  the  area  known  as  "Horseshoe  Park."  The  disturbed 
portions  of  these  streams  provide  abundant  sediment  sources  for  down- 
stream reaches  and  the  sediment  load  of  Fall  River  has  been  greatly 
increased. 

As  part  of  field  studies  begun  in  May,  1983,  a  sediment  sampling 
station  was  established  on  Fall  River  approximately  400  meters  down- 
stream of  the  alluvial  fan  terminus.  At  the  cross  section  selected,  the 
stream  is  relatively  straight,  has  near  vertical  banks,  uniform  bed 
topography  and  maintains  a  single  channel. 


Twenty-one  bedload  samples  were  collected  between  early  June  and 
mid-August  using  a  Helley-Smith  bedload  sampler  (Hel]ey  and  Smith, 
1971).  Measurements  were  taken  every  0.5  meter  across  the  channel 
giving  18  measurement  stations.  At  each  station,  the  sampler  was 
lowered  to  the  bed  for  a  fixed  time  period,  usually  30  seconds.  The 
collected  sediment  was  emptied  into  a  bucket  when  the  sampler  was 
approximately  two-thirds  full.  Thus  the  total  sample  actually 
represents  an  average  transport  rate  across  the  section.  Split  portions 
of  the  samples  were  sieved  by  the  U.S.  Bureau  of  Reclamation,  Denver, 
Colorado.  A  cumulative  frequency  curve  for  each  sediment  sample  was 
constructed  from  the  sieve  data.   Discharge  measurements  were  made  using 

a  Price  current  meter  and  wading  rod.   Flows  ranged  from  less  than  1.0 

3  3 

m  /sec  to  over  10.0  m  /sec  during  the  spring  snowmelt  period.   The 

duration  and  peak  of  runoff  in  1983  was  far  greater  than  the  long  term 
average. 

Hydraulic  data  obtained  during  sediment  sampling  or  stream  gaging 
included:  discharge,  width,  depth,  velocity  and  water  surface  slope. 
From  these  data,  values  of  cross-sectional  area,  unit  stream  power  and 
Froude  Number  were  calculated.  Unit  bedload  transport  rate  was 
calculated  from  the  total  weight  of  the  sample  collected.  Character- 
istics of  the  sediment  mixture  such  as  median  grain  size  (Dt-n)  and 
sediment  sorting  (a)  were  determined  graphically  from  the  cumulative 
frequency  curve  for  each  sample  (Folk,  1966).  The  data  used  in  the 
following  analyses  are  presented  in  Table  1. 

SEDIMENT  TRANSPORT  AND  HYDRAULIC  GEOMETRY  RELATIONS 

Sediment  characteristics  are  inextricably  tied  to  the  hydraulic 
properties  of  the  flow  in  problems  of  sediment  transport.   The  physical 


properties  of  the  sediment  determine  the  type  and  quantity  of  material 
to  he  transported.  The  hydraulic  properties  of  the  channel  and  the  flow 
define  the  transport  capacity  of  the  stream.  Simons  and  Richardson 
(1966)  suggest  that  in  the  simplest  sense  bedload  transport  rate  (i  ) 
can  be  represented  by: 

ib  =  f(Ys,Yf,T,v,D50,a,w,d,v,S)  (1) 

where  y  and  y  are  the  specific  weights  of  the  solid  and  fluid 
respectively,  T  is  the  water  temperature,  v  is  the  kinematic  viscosity 
of  water,  D  _n  is  a  characteristic  sediment  size,  <J  is  the  sediment 
sorting,  w  is  the  stream  width,  d  is  the  average  flow  depth,  v  is  the 
average  velocity,  and  S  is  the  water  surface  slope.  Many  of  these 
variables  are  not  independent  and  their  interactions  are  complex. 
Bedload  Transport 

Sediment  load  is  commonly  related  to  discharge  by  means  of 
"sediment  rating  curves."  Although  such  curves  offer  a  simple  predic- 
tive relationship,  sediment  transport  is  not  controlled  by  discharge. 
Bagnold  (1966)  first  introduced  the  concept  of  stream  power  per  unit  bed 
area  as  the  theoretically  correct  indicator  of  the  transport  ability  of 
the  stream.   Unit  stream  power  (u>)  as  defined  by  Bagnold  is: 

w  =  £gQS  =  yf  dSv  (2) 

w 

where  p  is  the  density  of  water,  g  is  the  gravitational  acceleration  and 

yf,  S,  d,  v,  and  w  are  as  defined  above.   Field  studies  by  Leopold  and 

Emmett  (1976)  have  verified  the  utility  of  this  approach. 

A  stepwise  linear  regression  procedure  was  used  to  develop  a 

multivariate  equation  to  express  bedload  transport  rate  as  a  function  of 

a  number  of  independent  variables.   In  the  stepwise  procedure,  variables 


are  added  to  the  regression  sequentially.  At  each  step,  the  regression 
is  re-examined  for  the  appropriateness  of  the  multivariate  model. 
Variables  with  non-significant  partial  F  values  are  removed  and  the 
model  is  refitted  using  the  remaining  variables.  Unit  stream  power  (uj)  , 
Froude  number  (Fr),  median  grain  size  (DrrJ  an^  sorting  (a)  were  used  as 
independent  variables.   The  resulting  regression  equation  is: 

log(ib)  =  -2.76  +  1.161og(w)  +  0.72(D5Q)  (3) 

where  i.  is  in  Kg/M/sec,  w  is  in  N/M/sec,  and  D^  is  in  millimeters. 
This  equation  has  a  coefficient  of  determination  of  0.67,  and  is 
represented  graphically  in  Figure  2.  Bedload  transport  rate  is  thus 
positively  correlated  with  stream  power  and  sediment  size.  The  positive 
correlation  between  transport  rate  and  unit  stream  power  is  as  expected 
but  the  positive  correlation  between  grain  size  and  transport  rate  is 
surprising.  Bedload  transport  rate  often  decreases  with  increasing 
grain  size  at  a  given  stream  power.  The  observed  relationship  suggests 
that  the  threshold  for  particle  motion  was  exceeded  for  most  of  the 
samples  taken. 
Hydraulic  Geometry 

The  hydraulic  geometry  of  river  channels  has  become  a  standard 
descriptive  tool  since  Leopold  and  Maddock  (1953)  first  introduced  the 
basic  relations.  They  demonstrated  that  changes  in  the  width,  depth  and 
velocity  of  flow  change  as  power  functions  of  discharge: 

w  =  aQb  d  =  cQf  v  =  kQm  (4) 

When  width,  depth  and  velocity  are  plotted  against  discharge  on  log-log 
paper  the  data  often  plot  as  a  straight  line.  However,  Richards  (1982) 
suggests  there  is  no  theoretical  basis  for  the  power  function  relation- 
ships and  there  may  be  significant  departures  from  the  straight  line 


model   where   changes   in   channel   roughness   occur  as   discharge   is 
increased . 

The  at-a-station  hydraulic  geometry  relations  for  the  Fall  River 
data  are  shown  in  Figures  3  and  4.  The  width-discharge  relation  is  not 
included  as  there  is  no  significant  change  in  width  with  change  in 
stage.   This  is  a  reflection  of  the  vertical  banks  of  the  cross  section. 

There  is  a  noticeable  non-linear  trend  in  both  plots  (i.e.,  Figures 
3  and  4)  and  the  data  are  best  fit  by  second  order  polynomials  rather 
than  straight  lines.  The  hydraulic  geometry  relationships  for  Fall 
River  are: 

log(v)  =  -0.60  +  l.ll(logQ)  -  0.45(logQ)2  (5) 

and 

log(d)  =  -0.37  +  O.lO(logQ)  +  0.40(logQ)2  (6) 

These  non-linear  trends  are  a  reflecton  of  the  influence  of  bedforms  on 
hydraulic  conditions.  In  a  sand  bedded  stream  the  bed  is  successively 
formed  into  ripples,  dunes,  plane  bed,  and  antidunes  as  stream  power  is 
increased  (Simons  and  Richardson,  1966).  The  development  of  bedforms 
gives  rise  to  a  particular  type  of  flow  resistance  known  as  "form  rough- 
ness." Dunes  offer  the  greatest  resistance  in  subcritical  (Fr  <  1.0) 
flows  of  the  type  observed  at  the  sample  site. 
Temporal  Variability 

Sediment  transport  in  mountain  streams  is  often  discontinuous 
reflecting  variability  in  sediment  supply  from  upstream  source  areas 
(Hayward,  1979).  The  trend  in  Fall  River  bedload  transport  rates  over 
the  summer  months  is  shown  in  Figure  5.  In  general,  bedload  transport 
rates  at  the  sample  site  follow  the  trend  of  the  snowmelt  hydrograph. 
However  several  departures  are  apparent.   Bedload  transport  rates  are 


low  in  the  early  part  of  June  even  though  discharges  are  moderately 
high.  In  early  July,  several  cases  of  high  transport  rates  for  moderate 
discharges  are  noted. 

Surveyed  channel  cross  sections  upstream  of  the  measuring  station 
show  episodic  scour  and  fill  (Figure  6).  Low  bedload  transport  rates 
correspond  with  periods  of  channel,  aggradation  as  the  upstream  reaches 
are  being  used  as  sediment  stores.  Subsequent  degradation  in  these 
reaches  results  in  high  measured  transport  rates  at  the  sampling 
station.  Thus  the  variability  in  bedload  transport  rates  is  a  function 
of  non-uniform  sediment  supply  from  upstream  reaches  rather  than  changes 
in  stage. 

CONCLUSIONS 

The  analyses  of  bedload  discharge  data  from  Fall  River,  Rocky 
Mountain  National  Park  can  be  summarized  as  follows: 

(1)  Stream  power  per  unit  bed  area  and  median  grain  size  are  good 

predictors  of  unit  bedload  transport  rate  as  indicated  by  the  good 

2 
correlation  (R  =  .67)  of  a  multiple  linear  regression. 

(2)  The  presence  of  bedforms  strongly  affects  the  hydraulic 
geometry  relationships  at  the  sampling  station. 

(3)  Periods  of  low  bedload  transport  rates  at  the  sampling  station 
correspond  to  episodes  of  channel  infilling  in  upstream  reaches  while 
high  transport  rates  correspond  to  episodes  of  scour. 
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Figure  2.   Bedload  transport  rate  (i  )  as  a  function  of  stream  power 
(u))  and  median  grain  size  (D   )  ; 
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At-a-station  hydraulic  geometry   relationship  between 
discharge  and  velocity. 
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Figure  4.      At-a-station  hydraulic   geometry   relationship  between 
discharge  and  depth. 
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As  the  Nation's  principal  conservation  agency,  the 
Department  of  the  Interior  has  responsibility  for  most  of 
our  nationally  owned  public  lands  and  natural  resources. 
This  includes  fostering  the  wisest  use  of  our  land  and 
water  resources,  protecting  our  fish  and  wildlife,  preserv- 
ing the  environment  and  cultural  value  of  our  national 
parks  and  historical  places,  and  providing  for  the  enjoy- 
ment of  life  through  outdoor  recreation.  The  Department 
assesses  our  energy  and  mineral  resources  and  works  to 
assure  that  their  development  is  in  the  best  interests  of 
all  our  people.  The  Department  also  has  a  major  responsi- 
bility for  American  Indian  reservation  communities  and  for 
people  who  live  in  island  territories  under  U.S.  admini- 
stration. 
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